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ABSTRACT 

The observational evidence for central black holes in globular clusters has been argued extensively, and their 
existence has important consequences for both the formation and evolution of the cluster. Most of the evidence 
comes from dynamical arguments, but the interpretation is difficult, given the short relaxation times and old ages 
of the clusters. One of the most robust signatures for the existence of a black hole is radio and/or X-ray emission. 
We observed three globular clusters, NGC6093 (M80), NGC6266 (M62), and NGC7078 (Ml 5), with the VLA in 
the A and C configuration with a 3-a noise of 36, 36 and 25 /iJy, respectively. We find no statistically-significant 
evidence for radio emission from the central region for any of the three clusters. NGC6266 shows a 2-cr detection. 
It is difficult to infer a mass from these upper limits due to uncertainty about the central gas density, accretion rate, 
and accretion model. 
Subject headings: globular clusters 



1. INTRODUCTION 

Although we do not understand how the nuclei of galaxies 
form or why they have black holes (BH) at their centers, the cor- 
relation between BH mass and bulge velocity dispersion does 
shed light on their formation and evolutionary histories (Geb- 
hardt et al. 2000a, 2000b: Ferrarese and Merritt 2000). A num- 
ber of different theories (e.g., Silk & Rees 1998; Haehnelt & 
Kauffmann 2000; Robertson et al. 2006) predict a BH mass 
bulge-velocity-dispersion relation, although they predict differ- 
ent slopes and intercepts for this relation. Exploration of the 
extreme ends of this relationship will help illuminate the un- 
derlying physical model, and in this paper we focus on the 
low mass end. Black holes at the low end of the relations, 
with masses between 100 and 10 6 M Q , are generally referred 
to as intermediate-mass black holes (IMBHs). There is signif- 
icant evidence that black hole masses from 10 5 -10 6 M Q exist 
from the work of Barth, Greene & Ho (2005) and Greene & Ho 
(2006). To go to yet smaller black hole masses, an extrapolation 
of the correlation between black hole mass and stellar velocity 
dispersion suggests studying stellar systems with velocity dis- 
persions of 10-20 km s" 1 . These dispersions are characteristic 
of globular clusters. Whether the existence of black holes in 
globular clusters could shed light on the formation and correla- 
tions of supermassive black holes is unknown, but clearly it is 
a possibility. Furthermore, the existence of massive black holes 
in clusters will have a significant effect on the cluster evolution. 
Thus, quantifying the demographics of black holes in clusters 
may be related to how supermassive black holes grow, and will 
definitely yield useful information about the evolution of clus- 
ters. 

Theoretical work suggests that we might expect IMBHs at 
the centers of steller systems (Ebisuzaki et al. 2001; Portegies 
Zwart & McMillian 2002; Miller & Hamilton 2002), although 
it appears to be difficult to make black holes more massive than 
100 M Q . Gurkan et al. (2004) suggest that IMBHs may be easy 
to form through runaway collisions with massive stars. Discov- 
eries of BHs in globular clusters have been claimed — Gl in 
M31 (Gebhardt, Rich & Ho 2002) and M15 (van der Marel 



et al. 2002; Gerssen et al. 2002). In fact, the M15 claim 
has been made for the past 30 years, starting with the result 
of Newell, da Costa & Norris (1976) and subsequently chal- 
lenged by Illingworth & King (1977). The basic issue is being 
able to distinguish a rise in the central mass-to-light ratio being 
due to either a black hole or the expected stellar remnants (neu- 
tron stars, massive white dwarfs and solar mass black holes). 
The most recent M15 result has been challenged by Baumgardt 
et al. (2003a). The result in Gl has also been challenged by 
Baumgardt et al. (2003b) but Gebhardt, Rich & Ho (2005) in- 
clude additional data and analysis that support the black hole 
interpretation. 

There has been two further observations which strongly sup- 
port the existence of a black hole in Gl. Trudolyubov & Pried- 
horsky (2004) measure X-rays from Gl using the Chandra Ob- 
servatory, centered to within 2" of the center of Gl. Subse- 
quently, Pooley & Rappaport (2006) suggest the X-ray emmis- 
sion is from accretion onto a black hole, and Maccarone & Ko- 
erding (2006) point out that if a black hole is present then a 30 
/iJy radio source may be expected. The most significant obser- 
vation comes from Ulvestad, Greene & Ho (2007) who find a 
28 /iJy (4.5er) emission centered on Gl. Other interpretations 
are a pulsar wind or a planetary nebula. The pulsar wind seems 
unlikely given the age of Gl and the point-like radio source (an 
old pulsar would have a large size). A planetary nebula would 
show optical emission lines which are not seen in the HST or 
Keck spectra of Gebhardt et al. (2003). 

Other studies of the existence of black holes in globular 
clusters have been less compelling. Colpi, Mapelli, & Pos- 
senti (2003) use indirect dynamical arguments to suggest a 
few hundred solar mass black hole in NGC 6752. McLaugh- 
lin et al. (2006) provide an estimate of black hole in 47Tuc 
of 900 ± 900 Mq . To date, there are no published upper lim- 
its of black hole masses that are significantly below that ex- 
pected from an extrapolation of the correlation between black 
hole mass and stellar velocity dispersion. 

While the dynamical arguements strongly support the black 
hole interpretation in at least Gl, the radio emission provides a 
clear and obvious result. Unfortunately, it is difficult to predict 
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the radio emission from a given black hole mass. The next step 
is to explore other globular clusters with a similar setup and 
deep exposures. 

2. TARGET SAMPLE AND RADIO FLUX DENSITY 
PREDICTIONS 

We selected three globular clusters for observation. First, 
using the stellar velocity dispersion at the center of Ml 5 
(NGC7078) suggests a revision of the mass of the BH at the 
cluster's center to 1000 M Q , higher than Maccarone's assumed 
value of 440 M Q , making M15 a promising candidate. Second, 
noting that Baumgardt et al. (2004) argue that highly centrally 
condensed globular clusters, as seen from their luminosity pro- 
files, are unlikely to harbor central IMBHs, we selected two 
globular clusters with large cores that are more likely to have 
central BHs. These clusters, NGC6093 and NGC6266, also 
have large central stellar velocity dispersions. 

In order to predict radio flux density, the first step is to use an 
expected black hole mass. The black hole masses can be esti- 
mated using an extrapolation of the correlations seen in galax- 
ies, namely either the black hole mass/velocity dispersion or the 
black hole mass/galaxy bulge luminosity relations. 

A precise prediction of the expected radio flux densities 
based the black hole mass is quite uncertain. Merloni et al. 
(2003) use radio flux densities, X-ray luminosities, and mea- 
sured black hole masses from both galactic and galaxian black 
holes to derive a fundamental plane for the three parameters. 
They argue that using any one parameter to predict another is 
quite uncertain. Unfortunately, X-ray luminosities do not exist 
for the three clusters studied here. Furthermore, the Merloni 
et al. study do not include any black holes with masses from 
10 to 10 6 , making any use of the fundamental plane suspect 
for the three globular clusters. Therefore, instead of directly 
using expected black hole mass and measured X-ray luminos- 
ity to predict the radio flux density, we simply use use the lo- 
cation between the 10 6 M Q black holes and the galactic black 
holes in the fundamental plane. In this region, the expected 
5 Ghz radio power ranges from 10 30 - 10 32 ergs/s. Indeed, for 
the Gl radio emission from Ulvestad, Greene & Ho (2007) cor- 
responds to 10 32 ergs/s, which is consistent with the measured 
black hole mass of 2 x 10 4 from Gebhardt, Rich, & Ho (2005). 
Thus, in order to predict expected flux densities, we adopt this 
range in radio power and use the known distances of the globu- 
lar clusters. A significant assumption in these estimates is that 
the physical conditions are similar; if, for example, the gas den- 
sity where much lower in globular clusters, the predicted radio 
power would be much less. 

Alternatively, Maccarone (2004) estimate expected radio 
emission based on the expected gas density and the correlation 
of Merloni et al. (2003). The gas density in the cluster come 
from the estimate of Freire et al. (2001) who use differences in 
column densities measured from pulsars in the front and back 
sides in the globular cluster 47 Tuc. While there is no reason 
to expect similar gas densities from cluster to cluster, it is the 
best measure we have of gas density in a cluster and therefore 
we adopt that value. Maccarone (2004) further assumes that 
the BH is accreting intra-cluster gas at 0.1 and 1% of the Bondi 
accretion rate. He assumes the BH mass to be 0. 1 % of the glob- 
ular cluster mass, which he estimates from the cluster's total lu- 
minosity and an assumed mass-to-light ratio, and computes the 
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expected 5 GHz flux density from the vicinity of the central BH 
for 15 globular clusters. Six of the globular clusters in Macca- 
roni's list lie north of the southern declination limit of the Very 
Large Arrajfl (VLA) and have an estimated 5 GHz flux density 
of 40 /iJy or greater (at 1% of the Bondi rate). We searched 
the VLA archive for observations of the centers of these clus- 
ters with noise levels low enough to have allowed a detection 
at Maccarone's predicted levels. No VLA archive data were 
found which had the required sensitivity. 

There have been two similar studies to the one presented 
here. Maccarone et al. (2005) provide upper limits for omega 
Cen using ATCA observations and for M15 from archival VLA 
observations. De Rijcke et al. (2006) provide upper limits for 
47 Tuc and NGC 6397 based on ATCA observations. 

Table 1 shows the computed flux densities for NGC6266, 
NGC7078, and NGC6093 at a frequency of 8.6 GHz. We as- 
sume a spectral index, a = -0.7 to provide predicted fluxes at 
8.6 GHz, the frequency of our observations. Our BH masses 
come from the BH/sigma correlation (shown in the second col- 
umn), and we also report those from Maccarone based on lumi- 
nosity (the first column). 

3. OBSERVATIONS 

Source positions, integration times on source, beam dimen- 
sions and position angles, and our 3cr limits are given in Table 
2. 

We used the position of the center of M15 determined by 
Noyola & Gebhardt (2006) using the optical surface bright- 
ness profile from the Hubble Space Telescope (HST), which are 
good to less than 1". We observed that position using the VLA 
for 7.5 hours on October 13, 2004 in the A configuration at 8.6 
GHz (A 3.5 cm), where the VLA has its maximum sensitivity. 
The resulting map has an rms noise level of 8.5 /iJy/beam and 
covers i=s 1 arcmin centered on the center of M15. We clearly 
see the source AC21 1 reported by Johnston, Kulkarni, and Goss 
(1991) about 1.5 arcsec northwest of the cluster center at with 
a peak flux density on our map of 144 //Jy. We do not detect 
the other known low mass X-ray binary, M15 X-2 (White & 
Angelini 2001), even at 1-sigma. The image also contains the 
planetary nebula K 648, for which we get: RA. 21:29:59.39, 
Dec. 12:10:26.46 (J2000). The measured flux density is 4.2 ± 
0.2 mJy and the deconvolved size is 1.5 x 0.7 arcsec. We do not 
see the pulsar PSR 2127+12A or any other point source near 
the center of M15 at or above a level of 25 /iJy /beam. 

The positions of the centers of NGC6093 and NGC6266 also 
come from Noyola & Gebhardt (2006) (with a similar accuracy 
of around 1"). These clusters were observed on August 11, 20, 
and 26, 2005 using the VLA in the C configuration at 8.6 GHz. 
Each map has an rms noise of 12 /iJy/beam. These maps cover 
an area of 5.4 arcmin centered on each cluster's center. No ra- 
dio source is seen at or above a level of 36 /iJy/beam near either 
cluster's center. In the case of NGC6093 a source is detected 82 
arcsec to the SE of the cluster center at R.A. 16:17:05.00, Dec. 
-22:59:47.3 (J2000) with a flux density of 0.32 ± 0.05 mJy. 
The distance from the cluster center makes it unlikely that this 
source is related to NGC6266. No source is seen convincingly 
in the image of NGC6266. 

Figure 1-3 plots the HST optical images (as used in Noyola 
& Gebhardt) with the VLA radio contours. None of the three 
clusters show a significant detection at 8Ghz at the cluster cen- 
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ter. However, NGC 6266 (M62) shows about a 2-a peak at the 
center. 




FIG. 1 . — HST optical image of M15 overlayed with the VLA contours of 
the central 10". Positive 1,2, and 3cr noise contours (8.5, 17, and 25.5 /ijy) 
are shown in green and negative are shown in red. The blue circle marks the 
center determined from Noyola & Gebhardt (2006), with a diameter of 0.5". 
North is up and East to the left. The radio source, AC21 1, is easily seen just 
north of the center. The other known X-ray source, M15 X-2, is not detected. 




FIG. 2.— HST optical image of NGC 6093 (M80) overlayed with the VLA 
contours of the central 10". Positive 1,2, and 3cr noise contours (12, 24, and 36 
/iJy)are shown in green and negative are shown in red. The blue circle marks 
the center determined from Noyola & Gebhardt (2006), with a diameter of 1". 
North is up and East to the left. 




FIG. 3.— HST optical image of NGC 6266 (M62) overlayed with the VLA 
contours of the central 10". Positive 1,2, and 3cr noise contours (12, 24, and 36 
/iJy)are shown in green and negative are shown in red. The blue circle marks 
the center determined from Noyola & Gebhardt (2006), with a diameter of 1". 
North is up and East to the left. There is about a 2-a positive 8Ghz signal at 
the center. 

4. DISCUSSION 

Failure to detect radio radiation at 8.6 GHz from the centers 
of three globular clusters does not prove that no globular clus- 
ters have IMBHs at their centers. Besides not having a black 
hole, other interpretations include 1) accretion by the BH could 
be episodic and we happened to observe the BHs in an "off- 
state", 2) the gas density could be much lower compared to 
galaxies, 3) the radiative efficiency may be lower than assumed 
(although the assumed efficiencies are already quite low), 4) or 
the accretion model may not be adequate in general. We would 
predict, using the relation of Merloni et al. (2003) or using 
standard accretion models and gas density estimates (as done in 
Maccarone 2004), that we should have detected radio radiation 
at 8.6 GHz if accretion is steady and the accretion rate times 
the Bondi rate is 10~ 4 x or higher. We would not have been able 
to detect the flux density predicted by a rate of 10~ 5 x or less. 
Ulvestad et al. (2007) estimate the fraction of the Bondi rate of 
just under 1% for Gl, but it is difficult to interpret due to the 
unknown radiative efficiency. For galactic black holes, the ra- 
diative efficiencies appear to vary greatly with some lower than 
10~ 5 (Lowenstein et al. 2001), although consistent with rates of 
around 10% of the Bondi rate. 

Models which predict 8.6 GHz flux densities from central 
BHs in globular clusters above about 25 /iJy/beam can be tested 
with the VLA currently. The EVLA should produce, for contin- 
uum observations, a sensitivity improvement of about a factor 
of 15, making 8.6 GHz flux densities above about 2^Jy/beam 
detectable. 
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Table 1 

8.6 GHz Flux Density Values 



Cluster 


M BH 


M BH 


Distance 


Flux Density 




Maccarone, M sun 


this paper, M sun 


kpc 




NGC6093 (M80) 




1600 


8 


2 x 10 J - 10 3 


NGC6266 (M62) 


450 


3000 


6 


3 x 10 3 -10 5 


NGC7078 (Ml 5) 


440 


1000 


10 


1 x 10 3 -10 5 



Table 2 
Observations 



Cluster 


RA 


DEC 


Integration 


Beam; Pos. Ang 


3 a Limit 




J2000 


J2000 


hours 


arcseconds; deg 


fi]y 


NGC6093 


16:17:05.00 


-22:59:47.3 


7.0 


3.9x2.3; -6 


36 


NGC6266 


17:01:12.96 


-30:06:46.2 


7.0 


4.7 x 2.2; -6 


36 


NGC7078 


21:29:58.35 


+12:10:01.5 


6.5 


0.2 x 0.2; -76 


25 



